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ABSTRACT 

During the recombination epoch, cosmic background photons couple not only 
to free electrons through Thompson scattering, but also to the neutral hydro- 
gen through Rayleigh scattering. This latter is ~ 2% effect for photons near 
the peak of the photon energy distribution at z = 800 and a ~ 0.2% effect at 
z = 1100. Including Rayleigh scattering in the calculation reduces Silk damping 
at fixed redshift, alters the position of the surface of last scattering and alters 
the propagation of acoustic waves. We estimate the amplitude of these effects. 
For the Microwave Anisotropy Probe (MAP), Rayleigh scattering increases the 
anisotropy spectrum by 0.1% at the most. For the highest frequencies of the 
Planck Surveyor, the effects of Rayleigh scattering are much more dramatic (de- 
creasing the anisotropy spectrum by 3% at ~ 550GHz and / ~ 1000). The 
relative difference between the spectra of low and high frequencies is imposed by 
an oscillation with a function of multipole / and the oscillation amplitude is up 
to 0.5% between 100 and 550 GHz. Rayleigh scattering also slows the decoupling 
between radiation and matter, but the effect is undetectably small. 



1. Introduction 

With the upcoming launch of the Microwave Anisotropy Probe (MAP)'^ in 2001 and the 
Planck Surveyor^ in 2007, the study of the cosmic microwave background (CMB) fluctuations 
is about to enter a new epoch. These satellites will characterize the statistical properties 
of the microwave background at better than the 1% level. This dramatic improvement in 
experimental capabilities demands that theorists attempt to achieve this level of accuracy 
in their calculations. Otherwise, systematic errors in the theoretical calculations will lead 
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to systematic errors in our interpretation. Over the past several years, several corrections 
to the CMB spectrum have been identified that systematically alter the predictions of a 
given model at the 1% level (Hu et al. 1995; Gnedin & Gnedin 1998; Seager, Sasselov & 
Scott 1999). Without these corrections, theoretical systematic biases could be as important 
as experimental error. This paper explores another possible source of systematic error: 
neglecting Rayleigh scattering in estimates of the CMB fiuctuation spectrum. 

In their seminal paper, Peebles & Yu (1970) note that Rayleigh scattering accounts 
for a few percent of the opacity near decoupling. To simplify their analysis, they ignore 
it and include only Thompson scattering in their calculation. Most subsequent analytical 
and numerical work has also included only Thompson scattering. This approximation makes 
the amplitude of cosmic microwave background fluctuations a simple function of frequency 
and greatly simphfies the equations. Hannestad (2001) considers the effect of Rayleigh 
scattering; however, since he does not do a frequency-dependent calculation, he significantly 
underestimates its importance. 

In this paper, we estimate the effects of Rayleigh scattering on microwave background 
fiuctuations. Our calculation begins with the analytical approach of Peebles (1980) and 
Hu & Sugiyama (1995): we treat the photon-baryon fiuid as tightly coupled through the 

surface of last scattering. While this approach is not exact, it enables us to estimate and 
identify the effects of Rayleigh scattering. In section 2, we discuss the propagation of acoustic 
waves in the baryon-photon ffuid. The additional opacity source reduces photon diffusion 
and also leads to a frequency dependence in the phase of the acoustic wave. In section 3, 
we identify four different observational effects. While some of these effects are important 
only at high frequencies, the suppression of Silk damping leads to observationally significant 
changes in the Rayleigh- Jeans regime. In section 4, we check our analytical estimates with 
a full numerical integration of the full Boltzmann equations using a modified version of 
CMBFAST. Section 5 concludes. In this paper, if not specified, the speed of fight is set to 



be 1. 



2. Acoustic Wave Propagation 



Rayleigh scattering adds an additional source of opacity to photon propagation: 




(1) 
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where Ug is the electron density, uh is the neutral hydrogen density, ax is the Thompson 
cross-section, and 

oo oo 

^eff = E - ^ E /y^r/ ^ (0•95ci?^)-^ for « (2) 

i=2 j=2 

where fij is the Lyman series oscillator strength, Uij is the Lyman series frequency, c is the 
speed of light and Ra is the Rydberg constant of hydrogen atoms (see Lang 1999, equation 
[1.306], and the oscillator strength in equation [2.118] and references therein). Looking at 
equation (1) and ignoring the frequency dependence, one sees that, roughly speaking, addi- 
tion of Rayleigh scattering corresponds to an increasing rig or D,f,h^ (where ^f, is the baryon 
density in units of critical density and the Hubble constant Hq — lOO/i km s"^ Mpc~^), 
which will be seen to be true at low frequencies (i.e., the increasing of anisotropy spectra at 
u < 150GHz) in §4. Other atomic and molecular processes are unimportant: the abundances 
of H2, HD and other molecules are low in the early universe (Stancil, Lepp & Dalgarno 1998) 
and neither He or H^ arc significant sources of opacity at the relevant frequencies. The frac- 
tional abundance of hydrogen in the n — 2 state peaks at 10~^^ just prior to decoupling 
(Stancil, Lepp & Dalgarno 1998), so Balmer hues are not a significant source of opacity. 

In this section, we explore the effects of a frequency dependent cross-section on acous- 
tic wave propagation. Assuming that the radiation is unpolarized, we begin by expand- 
ing the microwave background fluctuations into a smooth unperturbed piece f^^\y) = 
[exp{hpi//kBT) — (where hp is the Planck constant, /cb is the Boltzmann constant, u is 
the comoving frequency and T is the present temperature), and, in the conformal Newtonian 
gauge, a perturbation piece that obeys the linear collisional Boltzmann equation (Wilson & 
Silk 1981; Bond & Efstathiou 1984; Dodelson & Jubas 1995; Hu & Sugiyama 1995): 
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where = k • n, k is the comoving wavevector of the perturbation, n is the unit vector in 
the direction of propagation. A; = |k| is the comoving wavenumber, U}, is the baryon velocity, 
77 is the conformal time and z{ri) is the redshift at rj. The a is the scale factor and its value 
at the present time is oq. ^ and $ are the Newtonian potential and the perturbation to 
the intrinsic spatial curvature. The photon distribution function is expanded in Legendre 
polynomials: 

/(I) (i., fc, /.) = E //'^ k) {21 + l)Pi{ii). (4) 
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Ignoring changes in R{—3p^/4:pi„ where p-y and pb are the radiation density and baryon 
density respectively) and terms involving $ and ^, a moment expansion of equation (3) 
gives (e.g., Hu & Sugiyama 1996): 



ik r 



f^'\u,k) + 2n'>{u,k) 



--tkfi'\u,k) 



ao X[{1 + z{ri))u] 



ao X[(l + z(ri))u] 



(5) 



where we have truncated the expansion beyond the second moment and ignored the angular 
dependence of Compton and Rayleigh scattering. Overdots represent derivatives to the 
conformal time rj. The baryon velocity obeys the familiar mass and momentum conservation 
equations: 



5i){k) — ikub{k), 



(6) 



PbUb{k) = — 
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dlni/ 3 



iPbclkdb{k), 



where 5b is the baryon over-density and Cb is the baryon sound speed. This equation explicitly 
allows for the extra photon drag due to Compton and Rayleigh scattering, and, consequently, 
for the slightly later momentum decoupling and lower amplitude of the baryon fluctuations. 
But the heating on the baryon due to the extra photon drag of Rayleigh scattering can be 
neglected because of the small relevant ratio for comparing Rayleigh and Compton heating 
~ (me/mp)[(A"^/A(z/)) — 1] ~ 10^^, where mg/mp is the electron mass to proton mass ratio, 
(A"^/A(z/))(=A'^/A(?7), see equation 7) is the mean free path ratio of photons without Rayleigh 
scattering to photons with Rayleigh scattering, which can be estimated from Figure 1 later. 
We assume that the baryon sound velocity is small and ignore its pressure. 

We define a momentum weighted mean free path, which is equivalent to the Rosseland 
mean opacity in stellar atmospheres: 



1 
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Following Peebles (1980), we look for solutions proportional to exp(i J uj{k, rfj) . We expand 
the solutions in 5\ — (A[(l + z{rf))v\ — \{r])) 



fi'\u, k) - A{k)^^^^[l + h{k)8\ + o{8\)] exp 
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(8) 



- 5 - 



Substituting this back into equation (3), we find 



and (6): 

Ub{k) = -M{k)- . ^ ~ exp 

1 — luoykjRX 



iu!{k, r])dr] 



(10) 



We can insert equations (9) and (10) back in equation (5). Solving the equations for 
the frequency to second order in A yields, 

uj{k) — Cgk — i'jk'^X, (11) 



where 



' ■ (12) 
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+ 4{1 + R ) /5 
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The inclusion of Rayleigh scattering has altered the dispersion relation by reducing A due to 
the second term in equation (1). 

The temperature perturbation is defined by 5T/T = 0(z/, /i;,/i) = — f^^\u,k, jj) / ^^^^^^ 
and can be expanded in Legendre polynomials 0(i/, A;, /x) = Ylii=o ©/(^) k){2l + 1)P/ (/x). 

Thus, from equation (3), we have the evolution of temperature perturbation: 

e + ,^^6 + *) = -* + ^^^j^j^L_[e„ - e + ie.A(,) + (i4) 

If we ignore the angular dependence of Compton scattering, equation (14) has the 
solution: 

{Q + ^){u,k,ii,rio) = / [-(eo + * + /XK6)r(i/,7;,r;o) -$+^]e-"('^''''''°)e^'='^(''-''°)d7?, (15) 

Jo 

where rjo is the present epoch and r(z/, r/i, r/a) = ^ ^ a[(i+^(>;))^] ^^- 

Taking the multipole moments of equation (15) and setting Ub = 36i, we have for / > 2 
(Hu & Sugiyama 1995), 

k, r)o) ^ (Oo + ^){r).){-iyjiikAr),) + -^{-iy-^e,{r).)[lji_,{kAr).) 

rm 

-{I + l)j,+i(A;Ar;,)] + {-i)' / (^ - ^)ji{kAr,)dri, (16) 

J n» 
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where Ar^ = r)o — r), Ar)^ = r)o — r7*(j^), rjo is the present epoch and rj^liy) is the epoch 
at the last scattering. The fluctuations on the last scattering surface are obtained from 
(Go + = (00 + '$){ri^)D{k,u) and Qi{r]^) = Qi{ri^)D{k,u), where D{k,u) is the 

damping factor (see its definition in §3.1). The undamped WKB solution Qo{r]) and ©1(^7), as 
well as the potential "^{i]) and can be obtained from the Appendix in Hu & Sugiyama 

(1995). They are all assumed not to be affected by Rayleigh scattering here, which is 
confirmed to be reasonable by the consistency between our analytical and numerical results 
in §4. 

Integrating over all k modes of the perturbation, we have the power spectrum of the 
auto-correlation function: 



3. Observable Effects 

Rayleigh scattering produces several observable effects. We begin our discussion with 
effects that are detectable in the Rayleigh- Jeans part of the spectrum and then turn to effects 
that are detectable only in the Wien tail. 



3.1. Suppression of Silk Damping 

Our treatment of Silk damping follows the analytical approach developed by Hu & 
Sugiyama (1995). At a given wavenumber, the amplitude of Silk damping is the integrated 
photon diffusion length: 

kD{v)= f\kv')dv'- (18) 
Jo 

Since Rayleigh scattering reduces the photon-diffusion length, it reduces the amplitude of Silk 
damping. The damping factor at a given wavenumber is weighted by the photon visibility 
function: 

rVi) 

D{k,v) = I dr}V{v,r]) e'xjp 
Jo 

where the V{r]) is the photon visibility function 

V(z^, v) = -^{^, V, Vo) exp[-r (i/, rj)] , (20) 

where r is the photon mean free path (Figure 1). Figure 2 shows the damping factor as a 
function of k with and without Rayleigh scattering for a fiat vacuum-dominated universe 



(19) 
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with D,^ — 0.35, h — 0.65 and D,b — 0.05 (where is the mass density in units of the 
critical density). 

At low frequencies, the effect is small as photons decouple before Rayleigh scattering is 
important. At higher frequencies, the visibility function is shifted to lower redshifts where 
Silk damping is more important. This change in the damping factor produces a proportional 
change in the multipole spectrum at / ~ 6000A;. The right panels in Figures 3 and 4 show 
the amplitude of the changes produced by this effect. 



3.2. Frequency Dependent Surface of Last Scattering 

Since the photon cross-section is frequency-dependent, the location of the surface of last 
scattering will depend on frequency. Since photons emitted at lower redshift are emitted 
from regions with significant amounts of Silk damping, they contribute little to the observed 
temperature fluctuations. Because of this latter effect, the surface of last scattering is also 
k dependent: 

rvo „ 

{1 + z {k, I/)) = D{k, u)-^ / dr]—V{u, ri) exp 



k^ 



(21) 



3.3. Frequency Dependent Acoustic Wave Phase 

Including Rayleigh scattering, the sound horizon at the last scattering rf'^+^[r7*(/c, u)] — 
j^*^^''^^ Csdf] will be frequency dependent. The higher the frequency, the farther the sound 
horizon. It will shift the peak location by 5k,/k{v) = 1 — T^^~^^{r]^,)/rJ{r)^:), where the 
superscripts ^•^+^ and stand for the cases with and without Rayleigh scattering respectively. 
Thus, the peaks in anisotropy spectra will shift in the direction of decreasing /: 61/1 = Sk/k 
and ^ = '^Q^^i y. The higher the frequencies are, the more the peaks shift. 



3.4. Frequency Dependent Polarization Amplitude 

The amplitude of the microwave background polarization is proportional to the photon 
mean free path at the surface of last scattering (Zaldarriaga 1995). Since the photon 
mean free path depends upon frequency, the amplitude of the polarization is now frequency 
dependent: 
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Fig. 1. — The upper panel shows the ratio of the effective photon mean free path with 
Rayleigh scattering (A) to the photon mean free path without Rayleigh scattering (A"^). The 
lower panel shows that the photon visibility function y (i/, 77) shifts toward the present epoch 
with increasing frequencies. 
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Fig. 2. — Differences of damping factors: 5D / D{k, u) = D^'^^^{k, h')/D^{k, u) — 1. The su- 
perscripts and stand for the cases with and without Rayleigh scattering, respectively. 
At low frequencies, the increase of D^^~^^{k,i') comes mainly from the decrease of mean 
photon-difTusion length and the suppression of Silk damping; while at high frequencies, the 
decrease of D^^^'^{k, v) is mainly caused by the shift of visibility functions. 
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P {k, v) oc (A(i/, r])) = D{k, v)-^ / driV{v, rj) exp 
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In equation (22), we have weighted the photon mean free path by the visibihty function. 



(22) 



4. Numerical Integration 



In this section, we describe results of numerical simulations of the full Boltzmann equa- 
tions. We have revised the CMBFAST"^ code (Seljak & Zaldarriaga 1996) so that it can 
evolve the Boltzmann equations for frequency dependent scattering. The photon distribu- 
tion function is approximated by 50 different bins. The fluctuations at each frequency are 
evolved through the hierarchy of coupled equations (Wilson & Silk 1981; Bond & Efstathiou 
1984; Seljak & Zaldarriaga 1996): 

(23) 

A?i) , (24) 
, (25) 
-kjAi^i]j, for />2, (26) 

where kj{=—T) is the (now frequency dependent) differential optical depth and 
^Ti,j{=^i{^ji 1 / > 0) is the Z-th multipole moment of the temperature perturbation 

at frequency Vj. All other symbols and equations (e.g. = — = ^ and vi, — —iub) are 
as defined in Seljak & Zaldarriaga (1996). 

For the baryons, the momentum equation is now modified to include the frequency- 
dependent coupling to the photons: in Seljak & Zaldarriaga (1996) is modified: 

Vb = —Vb + clkSb + ^ XI /?'^j(3Atij ~ ^'') + (2"^) 

opb ^ 

where fj is the weight per bin. These modifications significantly slow CMBFAST, since 
it now has to evolve 50 times more variables. We have confirmed that this modified code 
reproduces standard results. 
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^CMBFAST is available at: littp://pliysics.nyu.edu/matiasz/CMBFAST/cmbfast.litnil. 
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Fig. 3. — Relative differences of anisotropy spectra with and without Rayleigh scattering: 
5Ci/Ci{v) — C^^^'^ {v) / Cf — 1 (solid curves). Basically, they keep the forms of their damping 
factor differences (see Fig. 2). The dotted curves represent the power spectra /(/ + l)Ci 
without Rayleigh scattering (in arbitrary unit). The oscillations of 5Ci/Ci show that the 
power spectra with Rayleigh scattering shift in the direction of decreasing /. The higher the 
frequencies are, the more the spectra shift (higher oscillation peak height). Note that the 
peak locations are not exactly same in the numerical results and in the analytical results. 
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Fig. 4. — Relative differences of anisotropy spectra witli and witliout Rayleigli scattering: 
6Ci/Ci{iy) = C^^^'^ {v) / Cf — 1. Tlie positive differences at low frequencies come from the 
decrease of mean photon-diffusion length. The differences between the numerical results and 
the analytical results are partly caused by slightly different peak locations (see Fig. 3). 
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In our numerical calculations, we use the recombination results of Ue and uh from 
RECFAST code to calculate the opacity to photon propagation (Seager, Sasselov & Scott 
1999). 

Our numerical calculations agree remarkably well with the analytical estimates of (§2). 
Figures 3 and 4 show their relative differences as a function of I and v. As seen from the 
figures, the analytical results are qualitatively consistent with the numerical results. At high 
multipoles and high frequencies, the anisotropy spectra in both results decrease because the 
visibility functions shift to lower redshifts where Silk damping is more important. The effect 
is up to 0.5% for v ~ 350GHz and 3% for v ~ 550GHz at / ~ 1000. At low frequencies 
and high multipoles, both the analytical and numerical results show that the anisotropy 
spectra with Rayleigh scattering are higher (~0.1%) than the spectrum without Rayleigh 
scattering, which is caused by the increase of damping factor D{k, v) at low frequencies 
(Fig. 2). At low I, that difference is too small to be physically significant because the cosmic 
variance uncertainty is comparably too large (~10% for / < 10). In Fig.3, comparing to the 
peaks of the anisotropy spectrum without Rayleigh scattering, all the oscillation peaks of the 
spectrum differences shift ~ 1/4 "period" in the direction of decreasing /, which supports 
our analyses about the phase shift caused by frequency-dependent sound horizons at the 
last scattering in §3.3. The oscillation amplitude of the spectrum difference is up to 0.5% at 
V ~ 550GHz and I ~ 1000. 

There are small difference between the numerical calculations and our analytical esti- 
mates. Note that the shift of the first peak is more evident in the analytical results. These 
differences may be due to the assumption that ©o, ©i and the potential are not affected 
by Rayleigh scattering, or they may be due to using equation (16) rather than the exact 
solution of equation (15). 

5. Discussion and Conclusion 

In this paper, we estimate the effects of Rayleigh scattering on the microwave back- 
ground. For MAP, Rayleigh scattering is a relatively minor correction. The largest effect 
is that it reduces the amount of Silk damping by ~ 0.1% at ~ lOOGHz and I ~ 1000. 
The systematic error is comparable to MAF's statistical errors. For the Planck Surveyor's 
highest frequencies, Rayleigh scattering can have significant effects on the shape and ampli- 
tude of the microwave background fiuctuation spectrum. Overall, it reduces amplitudes, C;, 
increasingly with increasing multipole number /, and, for observations at a single frequency, 
might be confounded with varying VLi,h'^. However the frequency dependence of Rayleigh 
scattering imparts a distinctive signature, so that the effects pointed at in this paper should 
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be easily corrected. At ~ 550 GHz, the anisotropy spectrum decreases as much as 3%. The 
peak shift between the spectra of 100 and 500GHz makes their relative difference oscillate 
as a function of /. The oscillation amplitude is as high as 0.5%, much larger than the sta- 
tistical uncertainties in the Planck Surveyor data. While the high-frequency channels with 
1/ > 353GHz are going to be used primarily for dust modeling, it is useful to note that the 
cosmic microwave background fluctuations will not be frequency independent. By looking at 
regions of low dust emission, the Planck Surveyor should be able to detect this effect. This 
detection would be an additional check of the recombination history of the universe, and a 
direct measurement of the formation of atoms in the early universe. 
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